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Observational constraints on the dark energy density evolution 
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We constrain tiie evolution of tlie dark energy density Q.d from Cosmic Microwave Background, 
Large Scale Structure, and Supernovae la measurements. While Supernovae la are most sensitive to 
the equation of state wq of dark energy today, the Cosmic Microwave Background and Large Scale 
Structure data best constrains the dark energy evolution at earlier times. For the parametrization 
used in our models, we find iuq < —0.8 and the dark energy fraction at very high redshift f2J^ < 0.03 
at 95 per cent confidence level. 
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Observations 
rious form of dark energy 
celerated expansion of our Universe. So far, the focus has 
been on the equation of state w = p/p of dark energy and 
in particular on its current value wq. However, if a dy- 
namical dark energy or quintessence arises from the time 
evolution of a scalar (cosmon) field one expects in gen- 
eral, that the equation of state changes as the Universe 
expands. Various parameterizations of w as a function 
of the scale factor a or redshift z have been investigated 
[II ElIIl in ■ Yet, it is rather the amount of dark 
energy fl^ than the equation of state (which is related to 
the derivative of fi^) that influenced our Universe in the 
past. In this spirit, suitably averaged quantities like flis 
[Tgl and Ogf have been used to describe the effects 
of dark energy at the time of last scattering and during 
structure formation. 

No parameterization of w or fi^ - no matter how com- 
plicated - will perfectly describe the evolution of dark en- 
ergy. Yet, some essential key features for a viable model 
seem to be the following: today, the amount of dark en- 
ergy should be ^ 70%. Going back in time, this value 
must have decreased considerably, as current constraints 
yield a fraction of dark energy at the time of last scatter- 
ing riis < 8% ^3 . Supernovae measurements tell us that 
this decrease must have occurred swiftly, as the slope of 
this decrease is reflected in wq < —0.7. 

Usually, observations at low redshift, such as Sne la 
measurements are combined with structure formation 
and CMB observations that are probing earlier epochs. 
In this paper, we will take a different point of view. Given 
the uncertainties in parameterizing f2rf or equivalently 
w, we look at high redshift and low redshift constraints 
separately. We use a particularly simple and direct pa- 
rameterization of the dark energy evolution |2(J|. The 
parameters are the amount of dark energy today fl^ and 
the amount of dark energy at early times 17^ to which 
it asymptotes for very large z. In terms of these, our 
parameterization is 
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FIG. 1: Constraints on the model parameters fij and 
wq. The left panel depicts the distribution from 
WMAP-hCBI-hVSA-hSDSS-HHST data, the right panel that 
of Sne la alone. The regions of 68% (95%) are enclosed by a 
white (black) line. The two data sets give almost orthogonal 
information; Sne la constrain wq but are less sensitive to 
while CMB and LSS are more sensitive to fij (and hence flis 
and Hsf) than to the precise value of wo- 
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The constant b is fully specified by the three parameters 
Wo, and fi^: 
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In terms of wq and b, the equation of state is w{a) — 
wq/{1 — 61na)^. For a comparison of our model to a cos- 
mological constant, see figures |21 and The advantage 
of a parameterization in terms of only three parameters 
is a minimal setting that accounts for information be- 
yond a Taylor expansion around z = (i.e. beyond the 
parameters and wq). Nevertheless, it embodies the 
most crucial potential new feature of a dynamical dark 
energy, namely the possibility of early dark energy. For 
supernovae, it seems more suitable than a continuation 
of the Taylor expansion which at the next step involves 
the derivative of the equation of state w'q — dw/dz \z=q- 
The present bounds on w'q allow a region of large \w'q\ 
for which the validity of a Taylor expansion is doubtful 
even for z = 1. 
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FIG. 2: Evolution of a typical dark energy model considered 
in this paper with 0.°^ = 0.7, wo = -0.999 and ^1"^ = 0.01 
depicted as dashed (red) line and with fijj — 0.03 depicted 
as dashed-dotted (blue) line. These two models correspond 
to values of f25 that are allowed within la and 2a respec- 
tively. For comparison, the solid (black) line shows the 
evolution for a cosmological constant with same 

Q.°a- The 

straight dashed-double dotted (orange) line at ild ~ 0.01 in- 
dicates the asymptotic limiting value of Qa for our model 
with n§ — 0.01. In contrast to the cosmological constant 
which contributes negligible at early times, our models have 
Qd > 0.01 (0.03) always and contribute a considerable effec- 
tive fraction Q^f ~ 0.06 (0.1) during structure formation. The 
dotted (green) line indicates the total equation of state of the 
Universe for our — 0.01 model. For Wtotai < —1/3, the 
Universe is accelerating. 



In addition to the dark energy parameters fl^, wq and 
ri^, we consider the matter and baryon densities today 
^m^^ ^^'^ ^6^^7 the Hubble parameter h, optical depth 
to the last scattering surface r and spectral index n. 

We compare the predictions of these models to the Sne 
la data of jlj] , as well as to the data from WMAP 
CBI 3, VSA [3, SDSS d and the Hubble parameter 
constraint of the Hubble Space Telescope 8] combined. 
For this, we employ the AnalyzeThis! 21j Monte Carlo 
package of CMBEASY We ran simulations both in 

terms of f2^, and wq as well as in terms of b, fi^ and wq, 
i.e we computed chains where the parameter fi^ is traded 
for the parameter b and vice versa. For a constraint on 
the parameter b from WMAP and Sne la, see Figure 01 

As far as the model parameters wq and are con- 
cerned, the main results are encapsulated in Figures ^ 
13 and As is seen from Figure ^ the information pro- 
vided by CMB plus LSS and Sne la are almost orthog- 
onal. While Sne la restrict wq considerably better than 
CMB and LSS combined, the sensitivity to H.^ is worse. 
From Sne la, we find wq < -0.78, while Q"^ < 0.029 from 
CMB, LSS and HST. 

The average dark energy fractions at last scattering 
and during structure formation, Clis and f2sf can substan- 
tially deviate from fl^ (as seen in Figures El and . As 
the underlying probability density is unknown, we cannot 
quote constraints on ilis and ilgf from counting the num- 
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FIG. 3: Constraints on the equation of state today. The dark 
(blue), medium (green) and light (red) shaded regions corre- 
spond to 1,2 and 3a confidence. The constraints from Sne 
la is displayed in the foreground (solid), while the less tight 
constraint from WMAP+CBI+VSA+SDSS+HST is depicted 
in the background. 



ber of models in the chain per bin. Instead, we choose to 
draw the distribution of the maximum likelihood model 
per bin. Figures and ^1 show this distribution. The 
somewhat tight 2 — a constraint on flis is a result of 
our parameterization: the amount of dark energy dur- 
ing structure formation Qgf always exceeds Clig. As the 
ISW effect in the CMB anisotropies as well as effects on 
the cold dark matter power spectrum restrict Clgf ~ 10%, 
this leads to fiis ~ 5% at 95% confidence level. 

Models of early dark energy typically lead to slightly 
lower values of h than A-CDM. The reason is the CMB: 
dark energy influences the CMB mostly due to projection 
effects and additional ISW contributions. To match ob- 
servations, any model must at least provide an acoustic 
scale I A that comes close to the one measured byWMAP. 
With fixed, the analytic expression for I a of 'iT] yields 



'^°^/o ^H"' i-no 



(4) 



where fir is the present energy fraction in radiation and 
w is the weighted average 
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The integral Q increases when w ^ — 1, i.e. for all other 
parameters fixed, the acoustic scale Ia becomes larger the 
more negative w. Conversely, as our models have w > 
— 1 by construction, we see that for all other parameters 
fixed, our models have a smaller acoustic scale compared 
to A-CDM. To counterbalance this, a somewhat smaller 
Hubble parameter h is preferred (see Figure 0), because 
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FIG. 4; Constraints on the parameter b from 
WMAP+CBI+VSA+SDSS+HST data (upper panel) 
and Sne la (lower panel). The dark (blue), medium (green) 
and light (red) shaded regions correspond to 1,2 and 3a 
confidence. 



I A depends on h to good approximation as |l 



Zyi cx 1 + /i 



1 + 0.4/1 ~ 



(6) 



where we have used the estimates w 0.7, ais w llOO"""^ 
and rij^g; h"^ « 4.4 x 10~^. Likewise, a sizeable early 
dark energy Oig can increase the acoustic scale accord- 
ing to cx 1/ y/l — fiis 01 . Both Cl\s and the somewhat 
smaller Hubble parameter counterbalance the effect of w 
in our models. 

Just as fiis expresses the main effect of early dark en- 
ergy on the CMB, the suitable average 



fist = [In atr. - Intteg.]" 
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with atr. = 1/3 encapsulates main effects of early dark 
energy on structure formation |l8j . By definition, one has 
f2sf(A) ~ 0.5% (i.e. non- vanishing) for a cosmological 
constant model and our choice of atr.- A sizeable flgf 
leads to a decrease in linear structure compared to A- 
CDM according to [3 
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FIG. 5: Supernovae compilation of Riess et. al. as data 
points with thin (brown) error bars. In addition we plo t the 
same data in binned fashion for illustration purposes 1241 . We 
plot the logarithm of the luminosity distance minus a fiducial 
model for which (IlHo = (l-l-z) In(l-l-z). The solid (black) line 
is for a cosmological constant model, the dotted (blue) line is 
for = 10"" and the dashed (red) line is for Q.'i = 10"\ AU 
models have wo = — 1 in common. 




FIG. 6: Likelihood distribution of the early dark energy pa- 
rameter inferred from WMAP-hCBI+VSA-|-SDSS-(-HST 
data. The dark (blue) , medium (green) and light (red) shaded 
regions correspond to 1,2 and 3(t confidence. 



Using Zeq. = 3500, we obtain for 3r2sf In Ze(j./5 ^ 1, the 
quick estimate 



(Js{D-E.) 
as (A) 



cx (l-6[af-at(A)]), 



(9) 



As we leave a free bias for our SDSS analysis, and have no 
prior on erg, one may worry about unphysically low values 
of (Tg for our models. It does turn out, however, that this 
is not necessarily the case. As seen in Figure |H1 our pre- 
dictions for cTg are compatible with observations, given 
the lower values of erg needed to explain non-linear struc- 
ture in early dark energy scenarios compared to ACDM 

m 

Cosmological probes have reached a level of accuracy 
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FIG. 7: Likelihood distribution for the hubble parameter h 
yielding h = 0.6±0.03. As explained in the text, slightly lower 
values of h than in standard ACDM models are preferred. The 
data used was WMAP+CBI+VSA+SDSS+HST. 
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FIG. 8: The linear rms power of matter fluctuations on 
8/i~^Mpc scales, ag. As ug is not a monte carlo parameter, 
we caution the reader that the underlying prior is not flat. 
With erg = 0.61 ± 0.08 the power is rather low compared to a 
standard A CDM model - a direct result of the suppression 
of growth due to the presence of dark energy during structure 
formation. 

that leads to stringent constraints on any model of our 
Universe. In the case of the parameterization of dark en- 
ergy in terms of fjj], fi^ and wq we used, the Supernovae 
la and CMB+LSS experiments provide almost orthogo- 
nal information. From the background evolution at z < 2 
that are probed by current Sne la observations, we find 
Wq < —0.78 at 95 % confidence level, in good agreement 
with prior investigations (e.g. |2,4 \2^). Likewise, struc- 
ture formation and the cosmic microwave background 
yield strong bounds on the abundance of dark energy 
at very high redshift. In our model, il^ < 0.029 at 95% 
confidence level. The value of fi^ is a factor of 2 — 3 
lower than the abundance of dark energy allowed during 
recombination Clis and the abundance during structure 
formation f2sf. This is both a feature and a shortcoming 
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FIG. 9: The average dark energy contribution during struc- 
ture formation Clis as defined in Equation ||7|l. As Hia is not 
a monte carlo parameter, but a derived quantity, we plot the 
maximum likelihood of all models per bin. In the model con- 
sidered, we see that a model with overall Likelihood less than 
Ax^ = 1 apart from the best fit model can be found only for 
Hi, < 3%. 




0.03 0.05 0.07 0.09 0.11 0.13 0.15 



FIG. 10: The average dark energy contribution during struc- 
ture formation Clsi as defined in Equation ||7|l. As Clsf is not 
a monte carlo parameter, but a derived quantity, we plot the 
maximum likelihood of all models per bin. Please note that 
Clsi ~ 0.5% by definition even for a cosmological constant 
Universe. 



of our simple parameterization. Curing it would require 
more parameters which at the present time are not nec- 
essary to explain current observations. As early dark 
energy models such as the one investigated in this paper 
make clear predictions for the CMB and in particular for 
structure formation, it will become possible to falsify or 
confirm such a scenario within the not so distant future. 
In this regard, the model investigated here holds more 
opportunities for detection than models of dark energy 
with negligible ild at higher redshifts. 
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